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Ratio of exit area to throat area 

Ratio of contour surface area to throat area 
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PERFECT BELL NOZZLE PARAMETRIC A PTIMIZATION CURVES 


INTRODUCTION 


Parametric nozzle performance data for bell nozzles are often required in the analysis and trade 
studies of rocket engines. In order to get accurate comparisons, it is very important to have a consistent 
set of data which covers the entire range of area ratios and nozzle lengths allowable. Although these 
data are available (in fragmented form) over a fairly wide range of area ratios, new high performance/ 
high pressure space engine applications require these data at area ratios well above those considered in 
the past. This document presents design data for un truncated nozzle expansion area ratios from 10 to 
6100 tor a specific heat ratio, (y) of 1.2 and includes curves which permit the optimization of nozzles 
for maximum thrust coefficient within a given length, surface area or area ratio. 


NOZZLE DESIGN 


The computer code documented in Reference 1 was used to perform the calculations for this 
report. The nozzle contours are constructed using the “method of characteristics” which is detailed in 
many texts such as Reference 2 and are of the wind tunnel or “perfect” nozzle type. These perfect 
bell nozzles are two-dimensional axisymmetric (i.e., all properties can be specified by two dimensions: 
axial distance downstream from the throat and radial distance from the centerline). They have the 
imposed constraint that the flow field velocity vectors at the nozzle exit (untruncated design area ratio, 
A d ) be of uniform magnitude and parallel to the centerline. Three thrust coefficients are calculated by 

the code for vacuum external conditions. They are the gross thrust coefficient, C fg , wliich is found by 

integrating the momentum and pressure force across the throat and then adding the thrust due to 
pressure on the nozzle wall; the net thrust coefficient, C fn , which is found by calculating the wall shear 

force (drag) caused by friction and subtracting it from the gross thrust; and the ideal thrust coefficient, 
Cfj, which is found by evaluating the one-dimensional isentropic flow equations at the given area ratio. 

A nozzle efficiency factor ETA n , is defined as Cf n /Cj- ( and is used to convert ideal thrust coefficients 

(and thereby one-dimensional isentropic specific impulse) to the actual two-dimensional axisymmetric 
case with friction. Since the wall friction is only applied to find the net thrust coefficient after the 
frictionless nozzle contour has been determined, it does not affect the contour or surface area. 

A description of the input required by the nozzle design code is shown in Table 1. The actual 
input values used to generate the data for the curves are shown in Table 2. Slight changes in the number 
of points on the centerline and exit radius were necessary to successfully execute the code for different 
ranges of the design area ratio, A[y. A corner expansion was chosen since it yields shorter nozzles than 

ones with a radius downstream of the throat for any given area ratio. Although the friction factor would 
depend on wall smoothness and the Reynolds number which is a function of combustion product density, 
nozzle size, and the viscosity of the combustion products, a value of 0.003 Fanning (0.012 Darcy- 
Weisbach) was chosen as representative. Large, high pressure nozzles would have a little lower and small 
low pressure nozzles a little higher friction factor. 
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The nozzle contours produced by the design procedures above are shown in Figure 1 where the 
radial distance from the centerline, Y, and the distance downstream from the throat, L, are normalized 
by the throat radius, R t . The untruncated design area ratio, Ap, represents the area ratio where uniform, 

parallel flow has been achieved for that contour. Since these nozzles are very long and derive very little 
thmst from th6 nearly parallel walls near their exit, the nozzles are nearly always truncated to a smaller 
area ratio as discussed later. The contours for Ap’s of 600 and larger are of course not shown to 

completion in this figure and the somewhat erratic values are the results of convergence tolerance as the 
code iterated to obtain the desired input values. It should be noted that only the Ap need be specified 

to define a nozzle contour and only two of the three variables Y/Rj (which isy^). L/R t and A^ need 

to be specified to fully define a truncated nozzle. To enable a comparison of the performance of the 
different contours, lines of constant net thrust coefficient, Cj n , (which includes friction) and wall surface 

area normalized to throat area, A s /A t , are superimposed. Although Figure 1 gives a good visualization 

of the nozzle contours and performance comparison it is sometimes difficult to accurately read data from 
it. 


Figures 2 through 8 give improved resolution and interpreting capability. Figures 9 and 10 
present the internal wall pressure normalized by chamber pressure ( i.e. , total pressure since the calcula- 
tions are for isentropic flow). This is the pressure that is integrated along the wall to find the wall thrust 1 

contribution. The important nozzle efficiency, which was discussed above in the nozzle design 

section is presented in Figures 11 through 13. Note that the efficiency has a maximum value for each 
area ratio. Since the one-dimensional, frictionless, ideal thiust coefficient is constant for a given area 
ratio, there is an Ap contour which produces a maximum thrust for a given area ratio. In other words, 

a set of nozzles exists which produce a maximum thrust for a given area ratio. This is covered further 
in the next section. 

NOZZLE OPTIMIZATION 


An Optimization Process is used to determine where to truncate the full nozzle contour to obtain 
maximum thrust ( i. e . , max Cf n ) within a given constraint such as nozzle area ratio (or diameter), surface 

area, or length. This optimization is best visualized by considering an enlarged section of Figure 1 shown 
in Figure A. Point “a” on the figure, which is the point where the net thrust coefficient, Cf n , line is 

tangent to a line of constant Y/R t (or area ratio, e). is the optimization representing the exit point of a 

nozzle contour yielding maximum thrust for a given area ratio or diameter. This is the optimization 
mentioned in the preceding section and as discussed there is also the point of maximum nozzle efficiency 
as seen in Figures 1 1 through 13. These points have been replotted on Figures 14 and 15 to enable 
nozzles of this set to be easily defined. Similarly, point b, the point at which the C'f n line is tangent to 

the constant surface area line, A g /A t , represents the optimization of maximum thrust for a given surface 

area. This set is given in Figures 16 and 17. Nozzles of maximum thrust for a given length are repre- 
sented by point C and curves for this set are shown in Figures 18 and 19. These minimum length 
optimized nozzle will be of a slightly different contour than the Rao minimum length nozzles [3| 
because the Rao nozzles do not have the added constraint which requires uniform parallel flow at the 
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untruncated design area ratio. The actual difference in performance at a given length will be very small 
between the two design methods, and the family of perfect bell nozzles has the advantage of including 
t data for the other optimized sets of minimum surface area and maximum efficiency for a given area 

ratio as well as any other nozzles up to the untruncated wind tunnel nozzle. Point “d” on Figure A 
represents the most thrust obtainable from any given nozzle contour, Ap. When the nozzle contour extends 

beyond this point, wall friction forces become greater than the pressure forces netting a negative thrust 
contribution which decreases the total thrust. This set of nozzles is not of practical interest, however, 
because the set of nozzles represented by point “a” are both shorter and have smaller diameters for the 
* same thrust. 


CONCLUSIONS 


Use of the optimization curves allows selection of nozzle contours with certain desired charac- 
teristics and comparison of nozzles of one optimum set to others. Minimum iength perfect Bell nozzles 
are similar in contour and performance to the minimum length nozzles designed by the Rao optimization 
method. 
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TABLE 1. BELL NOZZLE DESIGN PROGRAM INPUT DESCRIPTION 
PERFECT GAS MODEL OPTION 
GAMMA = 1.2 


Variable 

Default Value 

Description 

FCON 

1.0 

FCON = 1 .0 for two-dimensional 
axisymmetric flow 



FCON = 0.0 for two-dimensional plane flow 

R 

0.0 

Normalized radius of circle (R/R t ) for 
circular expansion 



R = 0.0 for corner expansion 

XNSL 

0.0 

Number of streamlines to be calculated 

CYLHT 

0.01 

Ratio of cylinder radius to throat radius 

XN1 

20.0 

Number of points on starting line 

XNE 

60.0 

Number of points on exit Mach line 

DLX 

0.1 : 

Maximum distance between adjacent Mach 
lines along the cylinder normalized by 
throat radius (AX/T t ) 

ATOL 

0.0005 

Iteration tolerance on design area ratio 

PITOL2 

0.005 

Maximum increment of tan 0 for corner or 
circular expansion 

EM 

1.005 

Starting line Mach number for uniform flow 

CF1 

0.003 

Incompressible coefficient of skin friction 
along nozzle contour 

DRAG 1 

0.011 

Tlmist loss at the throat 

TOLI 

0.00001 

Iteration tolerance for interior point 

PA 

0.0 

Ambient pressure ratio (Pa/Pc) 
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ABLE PERFECT BELL NOZZLE PARAMETRIC AND OPTIMIZATION ANALYSIS INPUT VALUES 
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Figure 3. L/D t versus Ap 
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Figure 8. THFTA versus A e /A t 
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Figure 14. Ap versus A e /A t for maximum ETA 





Figure 15. L/D t versus A e /A 
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Ae/ A t 

Figure 16. A D versus A e /Aj for minimum A s /A t nozzles. 
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Figure 1 7. L/D t versus A e /A t for minimum A /A t nozzles. 
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Figure 19. L/D t versus A c /A ( for minimum L/D { nozzles. 
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